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We highlight the emergence of metallic states in two-dimensional transition-metal-dichalcogenide 
nanostructures -nanoribbons, islands, and inversion domain boundaries- as a widespread and uni¬ 
versal phenomenon driven by the polar discontinuities occurring at their edges or boundaries. We 
show that such metallic states form one-dimensional wires of electrons or holes, with a free charge 
density that increases with the system size, up to complete screening of the polarization charge, and 
can also be controlled by the specihc edge or boundary configurations, e.g. through chemisorption of 
hydrogen or sulfur atoms at the edges. For triangular islands, local polar discontinuities occur even 
in the absence of a total dipole moment for the island and lead to an accumulation of free carriers 
close to the edges, providing a consistent explanation of previous experimental observations. To 
further stress the universal character of these mechanisms, we show that polar discontinuities give 
rise to metallic states also at inversion domain boundaries. These hndings underscore the potential 
of engineering transition-metal-dichalcogenide nanostructures for manifold applications in nano- and 
opto-electronics, spintronics, catalysis, and solar-energy harvesting. 


Transition metal dichalcogenides (TMDs) often crys¬ 
tallize in a layered structure so that two-dimensional (2D) 
mono- or few-layer sheets can be easily extracted from 
bulk materials, e.g. through exfoliation, similarly to what 
happens for graphene from graphite. The electronic and 
structural properties of such 2D TMDs -with chemical 
formula MX 2 - depend crucially on the transition metal 
atom M and on the chalcogen X (=S, Se, or Te) involved^. 
Recently, particular emphasis has been devoted to group- 
VI TMDs where M=Mo or W^’^. These 2D materials are 
semiconductors (with the only exception of the semimetal 
WTe 2 ) and offer the potential for extremely interesting 
technological applications^”^ in electronics, optoelectron¬ 
ics, (pseudo-)spintronics, photonics, and plasmonics. 

When designed into nanostructures like nanoribbons 
or triangular islands, group-VI TMDs become metal¬ 
lic at the edges, while remaining semiconducting in 
the bulk, owing to the appearance of edge states 
that cross the bulk energy gap. Such metallic edge 
states have been reported both theoretically^”^ ^ and 
experimentally^^’^^ along zigzag edges. Their pres¬ 
ence is associated with an enhanced photolumines¬ 
cence response^^ and with an extraordinary catalytic 
behavior^^, both in hydrodesulfurization^^’^^ and hydro¬ 
gen evolution reactions^^’^^. In addition, it has been the¬ 
oretically predicted that these metallic edge states are 
magnetic^”^^ with promising applications in spintronics 
and, more recently, that they give rise to unusual one¬ 
dimensional (ID) plasmonic excitations^^. Interestingly, 
there are both theoretical^^’^^ and experimental^^ evi¬ 
dences that similar metallic states appear also at inver¬ 
sion domain boundaries extending along a zigzag direc¬ 
tion. 

The existence of such metallic states at zigzag edges 
in group-VI TMD nanoribbons has been predicted to be 
a consequence of the non-zero bulk polarization of these 
materials^^’^^. Indeed, a polar discontinuity appears at 


the edges of nanoribbons with the appearance of localized 
polarization charges. These charges give rise to an elec¬ 
tric field that induces a charge reconstruction consisting 
in the accumulation of free carriers at the edge in order to 
screen the polarization charges and prevent a polar catas¬ 
trophe^^. This mechanism underlies similar insulator-to- 
metal transitions in SrTiOs/LaAlOs interfaces^^’^^ and 
their 2D counterparts^^’^^, and has important implica¬ 
tions in graphene-boron nitrid heterostructures^^. 

In this paper we address metallicity in group-VI TMD 
nanostructures in order to attain full control over their 
possible technological applications. With the help of ex¬ 
tensive first-principles density-functional-theory (DFT) 
calculations, we investigate how the charge density of 
free carriers varies with the system size and we explore 
the effects of specific edge/boundary configurations. In 
particular, we show that in some cases additional bound 
charges appear and affect the asymptotic value of the free 
charge density; as a paradigmatic example we investigate 
here the case of 2D TMD nanoribbons terminated with 
hydrogen or sulfur atoms, where metallic edge states are 
still presentWe show that the existence of a polar 
discontinuity is at the heart of the metallic states ap¬ 
pearing also in triangular islands and inversion domain 
boundaries. Last, such extensive investigation of the pe¬ 
culiar nature of these edge/boundary states paves the 
way for new exciting applications in solar-energy devices. 

Monolayer group-VI TMDs have an hexagonal lattice 
with three atoms (MX 2 ) in the unit cell. The atoms 
are arranged in three parallel monoatomic planes X-M- 
X, with transition-metal and chalcogen atoms forming 
two distinct 2D triangular sublattices (see Fig. 1). This 
gives each M atom a trigonal prismatic coordination and 
the overall point group is As already mentioned 

above, all group-VI TMDs are semiconductors (the only 
exception being the semimetal WTe 2 , for which the most 
stable structure is not the one described above), with a 
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FIG. 1: Panel (a): Top and lateral views of the crystal struc¬ 
ture of group-VI 2D TMDs (here M 0 S 2 ) with isosurface plots 
of the two kinds of valence Wannier functions: (i) a p-like 
Wannier function centered on a chalcogen atom (yellow), hy¬ 
bridized with transition-metal d-orbitals; and (ii) a Wan¬ 
nier function centered in the middle of the hexagonal cell, 
with contributions from the d-orbitals of the transition metal 
atoms (black) surrounding it. Panel (b): Schematic represen¬ 
tation of the crystal honeycomb structure where in addition 
to chalcogen (yellow) and transition metal (black) atoms also 
the centers of the Wannier function (ii) above are reported 
(small red circles). A non-primitive rectangular unit cell is 
emphasized with gray shading. When repeated periodically 
along the p-axis, this structure give rise to a AV = 4 zigzag 
nanoribbon, while when repeated along x it leads to a = 4 
armchair nanoribbon. 


direct band gap at the two inequivalent corners of the 
Brillouin zone. The bulk polarization^^’^^ of 2D TMDs 
can be conveniently computed by mapping the ground 
state of the system into a set of maximally-localized Wan¬ 
nier functions^^. Then, the electronic contribution to the 
polarization is a sum over point-like charges located at 
the Wannier centers {r)j and the total polarization reads 

/ N iVei/2 \ 

\a=l j=l J 

Here and Tq, are the charges and positions of the N 


ions in the unit cell, A^ei is the number of electrons, and 
E is the area of a unit cell. In Eq. (1) is a generic 
Bravais lattice vector, since in the modern theory of 
polarization^^the formal polarization is not a single 
vector (like in classical electrostatics) but rather a lattice 
of vectors. If we consider as valence electrons only the 
outermost p-electrons of the chalcogens and the s- and 
d-electrons of the transition metal, we then have that 
Zx = 4, Zm = 6, and A^ei = 14. In the following, max¬ 
imally localized Wannier functions corresponding to the 
valence bands of 2D TMDs have been computed using the 
WannierQO code^^ starting from the DFT ground state 
obtained with Quantum-ESPRESSO^^ (see Methods for 
more details). Owing to their isovalence and isostruc- 
tural properties, the valence Wannier functions have the 
same character for all group-VI 2D TMDs (in the case of 
WTe 2 we consider a metastable semiconducting structure 
identical to that of the other group-VI TMDs). In par¬ 
ticular, we find three p-like Wannier functions centered 
on each chalcogen, hybridized with transition-metal d- 
orbitals, together with a Wannier function centered in 
the middle of the hexagonal cell, with contributions from 
the d-orbitals of the atoms surrounding the cell^^ (see 
Eig. 1(a)). With this set of valence Wannier functions 
and their centers it is straightforward to compute the 
bulk formal polarization according to Eq. (1). Eor conve¬ 
nience we consider a non-primitive rectangular unit cell, 
as shown in Eig. 1, with the long side along an armchair 
direction (x-axis) and the short side along a zigzag direc¬ 
tion (^-axis). The bulk formal polarization is then given 

by 


P = 


2e . 

— ■ 

3a 


( 2 ) 


where a is the lattice constant and we set = 0 in 
Eq. (1), thus selecting a specific representative element 
of the polarization lattice associated with the choice of 
the unit cell in Eig. 1^6,37 stress that P is quan¬ 

tized and points along an armchair direction, in agree¬ 
ment with the constraints imposed by the three-fold ro¬ 
tation symmetry^^’^^. 

In a nanoribbon the periodicity of the bulk crystal is 
preserved along a given direction, which defines the di¬ 
rection of the nanoribbon, while it is broken in the or¬ 
thogonal direction where the edges create a discontinuity. 
Since vacuum can be considered as an insulator with zero 
polarization, we have that across the edge of a group-VI 
TMD nanoribbon a polar discontinuity arises. Polariza¬ 
tion charges will thus appear with a linear charge density 
given by 


Ap = P-n, (3) 

where n is a unit vector orthogonal to the edge. By com¬ 
bining Eq. (2) and (3), we see that when a nanoribbon 
extends along a zigzag direction, the polarization is or¬ 
thogonal to the edge, so that the polarization charge is 
maximal and it is given by Ap = ±2e/(3a), where the 
plus sign refers to the Mo-edge while the minus sign to 
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FIG. 2: Panel (a): Band structures of zigzag M 0 S 2 nanorib¬ 
bons with different widths: = 1,2, and 14. The dashed 

line highlights the position of the Fermi energy. Panel (b): 
Macroscopic average of the electrostatic potential energy de¬ 
fined in Eq. (4) as a function of the coordinate across the 
ribbon. Results are reported for nanoribbon widths = 
1,2,3,4,6,8,10,12,14. 


the S-edge. For armchair nanoribbons instead the polar¬ 
ization is parallel to the edge and the polarization charge 
vanishes (Ap = 0). This difference in the relative ori¬ 
entation between the edges and the bulk polarization 
is at the heart of the different electronic properties re¬ 
ported in the literature for zigzag and armchair TMD 
nanoribbons^’^’^“^^. Indeed, as we will show below, it is 
the electric field associated with the polarization charges 
that triggers an insulator-to-metal transition and drives 
the appearance of metallic edge states in zigzag (but not 
in armchair) nanoribbons. 

Let us first consider zigzag nanoribbons, obtained 
by the repetition of the rectangular unit cell shown in 
Fig. 1(b) along the x-direction times, while being 
periodic along the ^-direction. In Fig. 2(a) we report 
the band structure of M 0 S 2 nanoribbons with different 
widths, corresponding to = 1,2, and 14, respectively. 
Qualitatively similar results can be obtained for all other 
2D TMDs. For very narrow nanoribbons (A^^(;=l) the 
system is semiconducting, although the energy band gap 


is significantly reduced with respect to the bulk owing 
to the presence of mid-gap states that are localized at 
the edges. For larger ribbon widths, the gap closes and 
the system becomes metallic with a charge transfer from 
the top valence band to the bottom conduction bands. 
To understand this insulator-to-metal transition, we also 
show in Fig. 2(b) the macroscopic average^^ v{x) of the 
electrostatic potential energy v{x^y^z)^ defined as 


2 nX-\--s/^a/A na noo 

v(x) = — / / / v{x — x'z)dxdydz. (4) 

^ Jx — \/3aj4: Jo J —00 


The presence of polarization charges at the edges is con¬ 
firmed by the change in slope of v{x) across each edge. 
Indeed, although this slope is not related to the actual 
electric field acting on the electrons owing to the integral 
over an infinite range in the definition (4), Gauss theo¬ 
rem still applies and the change in slope of v{x) gives an 
important information on the presence of a finite amount 
of total charge at the nanoribbon boundaries and thus on 
the existence of polarization charges. These polarization 
charges induce an electric field that in turn gives rise to 
a potential energy difference Av between the two edges 
of the nanoribbon. This means that valence states on 
one edge are pushed upwards in energy, while conduction 
states at the opposite edge are pushed downwards, with 
an overall reduction of the global energy gap of the sys¬ 
tem. As the width of the nanoribbon increases, Av grows 
logarithmically^^, with a progressive decrease of the gap. 
Above a critical width, the gap vanishes and the system 
becomes metallic, thus explaining the evolution of the 
band structure provided in Fig. 2(a). This insulator-to- 
metal transition as a function of the nanoribbon width is 
associated with a charge reconstruction: electrons from 
the top valence bands are transferred to the bottom con¬ 
duction bands with the creation of electron and hole 
pockets of free carriers. In Fig. 3(a) we show the spatial 
profile of the charge density of free carriers for a = 6 
M 0 S 2 nanoribbon, as obtained by integrating the local 
density of states associated with the pockets of free elec¬ 
trons and holes. We note that these free carriers form 
ID metallic wires extending along the nanoribbon edges 
and present opposite character (electron or hole) at op¬ 
posite edges: electrons on the Mo-edge and holes on the 
S-edge. In the following we define as Ap the total free 
charge (per unit length) present in such ID wires, spec¬ 
ifying wherever necessary if we refer to the chalcogen or 
transition-metal edge. 

Further evidence of the electrostatic origin of metallic 
edge states is given in Fig. 4, where we show the band 
structures of M 0 S 2 zigzag nanoribbons with different 
widths obtained using a tight-binding model. The tight- 
binding parameters have been obtained for a bulk mono- 
layer M 0 S 2 by mapping both the valence and the con¬ 
duction bands into a set of atom-centered Wannier func¬ 
tions corresponding to Mo d-orbitals and S p-orbitals^^. 
Although edge states are still present owing to the under¬ 
coordination of edge atoms, the system remains always 
semiconducting with a gap that does not change with the 
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FIG. 3: Real-space distribution of the charge density of free 
electrons (blue) and holes (red) for a = 6 zigzag M 0 S 2 
nanoribbon with (a) bare and (b) hydrogen-terminated edges. 


nanoribbon width. As a consequence, in this case there 
is no accumulation of free carriers at the edges. Since 
electrostatic effects are not included within this tight- 
binding description, this means that, although polariza¬ 
tion charges and their corresponding electric field are not 
necessary to account for the existence of edge states in 
TMD nanoribbons, they are needed to drive metallicity. 
In addition, comparing in more detail the band struc¬ 
tures reported in Figs. 2 and 4 we observe that the elec¬ 
tric field leads also to a relative shift between the energy 
bands. This effect is particularly relevant on the Mo-edge 
where the electrons that, according to the bulk calcula¬ 
tions, should occupy the Wannier function centered out¬ 
side the ribbon (see also Fig. 1(b)) feel the attraction of 
the positive polarization charges and are thus transferred 
on the Mo-atom. This brings a band, which according 
to the tight-binding model would be unoccupied, below 
the Fermi energy and correspondingly a tight-binding va¬ 
lence band is pushed above the Fermi energy. The charge 
reorganization associated with the change in the band 
occupation can be effectively described in terms of an 
additional -f2e bound charge at the original position of 
the Wannier center and a — 2e charge on the external Mo 
atom, so that the overall charge at the edge remains the 
same. 

It is important to mention that the free charges asso¬ 
ciated with the edge states partially screen the polariza¬ 
tion charges, so that the total amount of charge at each 
edge decreases when the system becomes metallic. In¬ 
deed, in Fig. 2(b) we observe that the slope of v{x) drops 
when the nanoribbon from semiconducting {N^j = 1) be¬ 
comes metallic {N^ =2). In addition, the slope further 



FIG. 4: Band structures of zigzag M 0 S 2 nanoribbons with dif¬ 
ferent widths (Aw = 1,2, and 14). Differently from Fig. 2(a), 
in this case the energy bands are obtained from a tight- 
binding model that does not include electrostatic effects. 


decreases as the nanoribbon width increases, suggesting 
that the charge density of free carriers Ap grows with 
Nyj. This is due to the fact that, as the nanoribbon 
width increases, the electric field drives more and more 
free charges to the edges, until (asymptotically) these 
perfectly screen the polarization charges and the elec¬ 
tric field vanishes. To further confirm this picture, in 
Fig. 5 we show |Af| at the chalcogen edge as a function 
of the nanoribbon width for all group-VI TMD. In agree¬ 
ment with the above analysis on the slope of F(x), it can 
be seen that the free charge increases with the nanorib¬ 
bon width and asymptotically it equals the polarization 
charge, jApj = |Ap| = 2e/(3a). Qualitatively similar be¬ 
havior can be observed for all group-VI TMDs, with only 
minor quantitative differences in the rate at which the 
asymptotic limit is reached. Differences emerge instead 
when considering the free charge on the transition metal 
edge. Indeed, while in all Mo-based TMD nanoribbons 
the free charge on the Mo-edge is identical (in magnitude) 
to the one on the chalcogen edge, for W-based nanorib¬ 
bons we have observed that the free charge on the W-edge 
exceeds by 2e/a the free charge on the chalcogen edge for 
all widths. This is due to the fact that also the band as¬ 
sociated with the extra electrons on the external W-atom 
becomes metallic in WX 2 nanoribbons. 

In Fig. 6 we summarize what happens for armchair 
nanoribbons that result from the A^^-fold repetition of 
the rectangular unit cell shown in Fig. 1 along the y- 
direction, while preserving the bulk periodicity along the 
x-direction. In agreement with previous literature^’^~^^ 
we find that such nanoribbons remain semiconducting at 
all widths, although the energy gap is reduced with re¬ 
spect to the bulk owing to the presence of mid-gap states 
associated with the undercoordination of the edge atoms. 
Indeed, in this case an insulator-to-metal transition does 
not occur owing to the absence of polarization charges 
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FIG. 5: Charge density of free carriers jAp | at the chalcogen- 
edge as a function of the nanoribbon width for: (a) Mo-based 
and (b) W-based group-VI transition metal dichalcogenides. 
Different symbols identify different chalcogen atoms: S (black 
circles), Se (blue squares), and Te (red triangles). In all 
cases the free charge density asymptotically reaches the value 
2e/(3a) (dashed line), corresponding to perfect screening of 
the polarization charge Ap. Lines are hts according to Ref. 24. 
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FIG. 6: (a) Energy band structure and (b) profile of the 
macroscopic average of the electrostatic potential energy for 
a Vw = 8 M 0 S 2 armchair nanoribbon. 

at the edges. This is most evident in Fig. 6 where we 
see how the macroscopic average of the electrostatic po¬ 
tential v{x) remains flat inside the nanoribbon, meaning 
that no net charge is present at the edge. 

Up to this point we have assumed perfect zigzag or 
armchair ribbons where the edge atoms remain underco¬ 
ordinated. In the following we would like to address ter¬ 
mination effects, when additional atoms are chemically 
bonded at the zigzag edges to restore bulk coordination. 
For definiteness we focus on M 0 S 2 and consider hydro¬ 
gen or sulfur as saturating the edge dangling bonds, since 


these are species that are typically present during growth. 
In the case of hydrogen, we assume to have two H atoms 
bonded to the external molybdenum on the Mo-edge in 
order to recover bulk six-fold coordination, while at the 
S-edge each sulfur atom binds to a single H atom. We 
consider separately each edge since these can be assumed 
decoupled for sufficiently large nanoribbon widths. At 
the Mo-edge, in order to form a bond with the hydro¬ 
gen atoms two electrons are needed. This means that 
the total amount of bound charge A^° at the Mo-edge 
is not simply given by the polarization charge Ap°, but 
there must be an additional —2e contribution per unit 
length: A^° = Ap° — 2e/a = —4e/(3a). Similarly, on 
the S-edge the S-H bond is partially ionic, with sulfur 
removing an electron from the hydrogen. As a conse¬ 
quence, on this edge there is an additional +2e bound 
charge per unit length due to the hydrogen ions, and 
the total bound charge on the S-edge is thus given by 
A^ = Ap + 2e/a = 4e/(3a). Hence, we notice that with 
respect to pristine edges hydrogen termination leads to 
a change in the sign of the bound charge on both edges. 
For narrow nanoribbons, these bound charges are not 
compensated by free charges and create an electric field 
which points in the opposite direction with respect to 
pristine edges. In Fig. 7(b) we show the macroscopic 
average of the electrostatic potential energy v{x) (blue 
line), confirming that the slope is inverted with respect to 
the case of pristine edges (black line), signaling a change 
in the sign of the bound charges in agreement with the 
discussion above. The ensuing electric field will drive an 
insulator-to-metal transition as it was the case for pris¬ 
tine edges, and for sufficiently large ribbons the system 
becomes metallic, as shown in Fig. 7(a) where we plot 
the band structure for a M 0 S 2 zigzag nanoribbon with 
hydrogen-terminated edges. The only difference with re¬ 
spect to the pristine case is that now there are free holes 
on the Mo-edge and free electrons on the S-edge. This 
is reported in Fig. 3(b) where we plot the charge den¬ 
sity of free carriers for hydrogenated edges, and show 
that electrons and holes are exchanged with respect to 
Fig. 3(a). Asymptotically, the free charge has to screen 
the bound charge in order to prevent a divergency in the 
electrostatic energy. In Fig. 7(c) we show the free charge 
density |Af| as a function of the nanoribbon width con¬ 
firming that asymptotically there is perfect screening of 
the bound charge: |Af| ^ = l^bl ~ 4e/(3a). In 

summary, we have shown that the presence of hydro¬ 
gen on both edges does not prevent the occurrence of an 
insulator-to-metal transition and of a finite electric field 
inside the nanoribbon, although its sign is reversed with 
respect to pristine edges and consequently also the char¬ 
acter of free carriers along the ID wires at the edges is 
exchanged. 

Before discussing the case of S-termination we con¬ 
sider the case in which hydrogen atoms are added only 
at the Mo- or at the S-edge. According to the analy¬ 
sis above, in these cases we have that the bound charges 
have the same sign on both edges but different values. 
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FIG. 7: Panel (a): Energy bands for M 0 S 2 nanoribbons (iVw = 6) with different edge terminations: bare edges (black), 
hydrogen-terminated edges (blue), hydrogen-terminated S-edge (green), and hydrogen-terminated Mo-edge (red). Panel (b): 
Macroscopic average of the electrostatic potential energy defined in Eq. (4) as a function of the coordinate across the ribbons 
introduced in panel (a). Panel (c): Charge density of free carriers as a function of the nanoribbon width for different edge 
terminations. Large empty circles represent the charge density on the S-edge, |A||, while small filled circles refer to the Mo-edge, 
|A^°|. Depending on the edge termination, two different asymptotic limits are reached, at 2e/(3a) or 4e/(3a) (dashed lines). 
Lines are hts according to Ref. 24. The color coding in panels (b) and (c) is the same as that in panel (a). 


For instance, when hydrogen is present only at the Mo- 
edge we have that the bound charge is always nega¬ 
tive and given by A^° = —4e/(3a) on the Mo-edge and 
by A® = -2e/(3a) on the S-edge. Analogously, when 
only the S-edge is terminated with hydrogen, the bound 
charge is always positive and is A^° = 2e/(3a) on the Mo- 
edge and A^ = 4e/(3a) on the S-edge. In order to screen 
the electric fields associated with these bound charges, 
a charge reconstruction occurs and free carriers appear 
with the same character on both edges. The system is 
thus metallic, as shown in Fig. 7(a) where we report the 
band structures both in the case of hydrogen only on the 
Mo- or on the S-edge. In Fig. 7(c) we show the charge 
density of free carriers as a function of the nanoribbon 
width and we note that for sufficiently large widths the 
charge density of free carriers is constant and identical 
to the bound charge density and it is thus different on 
the two edges. This means in particular that at such 
nanoribbon widths the total net charge at each edge is 
zero and there is not an overall electric field. Indeed, in 
Fig. 7(b) the slope of v{x) stays flat across the ribbon 
when hydrogen is present only on the S-edge (green line) 
or on the Mo-edge (red line). Thus, in the hybrid situa¬ 
tion in which only one edge is terminated with hydrogen 
we no longer have an overall electric field across the rib¬ 
bon but we still have metallic edge states. Only for very 
small nanoribbon widths hybridization between the two 
edges gives rise to a charge transfer with a slight devia¬ 
tion of the free charge density from its asymptotic limit 
and thus with the appearance of a small electric field. 

We now turn our attention to the case when additional 
sulfur atoms are present on the Mo-edge while the S-edge 
remains pristine. This situation is extremely relevant for 


experiments since the bare Mo-edge is not stable over 
a wide range of sulfur chemical potentials^’^^, while the 
bare S-edge is favored with respect to edge hydrogenation 
or desulfurization^. In order to restore the bulk six-fold 
coordination of Mo atoms and stabilize the Mo-edge, it is 
energetically favorable to bind additional S atoms. Two 
plausible configurations have been identified^^’^^, which 
differ in the sulfur content at the edge, i.e. in the number 
of S atoms adsorbed per unit length. One is the S dimer 
configuration that results from the adsorption of two ad¬ 
ditional S atoms per unit length that bind to the same 
Mo atom, similarly to what happens in the bulk. Another 
is the S monomer configuration that corresponds to the 
adsorption of a single S atom per unit length, which lies 
in the Mo-plane and binds to two neighboring external 
Mo atoms. The sulfur content at the edge with respect to 
the bulk is thus 100% for S dimers, 50% for S monomers, 
and 0% for bare edges. 

For S monomers, since the S atom has a larger electron 
affinity than Mo, its presence modifies the edge charge 
reorganization described above for bare edges. In this 
case the two electrons that would be associated with a 
Wannier function centered outside the ribbon are taken 
by the S monomer rather than by the outer Mo atom. 
This does not change the total amount of bound charge, 
which is then simply given by the polarization charge, 
A^° = Ap°, as for bare edges. As a consequence, there 
is still an electric field that drives an insulator-to-metal 
transition and asymptotically the charge density of free 
carriers will be the same as for bare edges. This is con¬ 
firmed by the results of Fig. 8 where it is shown that: 
(i) nanoribbons with S-monomer terminated Mo-edges 
are still metallic; (ii) from the finite slope of v{x) one 


























7 





Width [A] 

FIG. 8: Panel (a): Energy band structures for M 0 S 2 nanorib¬ 
bons (A^w = 6) with different sulfur content on the Mo-edge: 
bare edge (black), S monomers (blue) and S dimers (green). 
Panel (b): Macroscopic average of the electrostatic potential 
energy defined in Eq. (4) as a function of the coordinate across 
the ribbons introduced in panel (a). Panel (c): Charge den¬ 
sity of free carriers on the S-edge, |A||, as a function of the 
nanoribbon width for different Mo-edge terminations. Lines 
are fits according to Ref. 24. The color coding in panels (b) 
and (c) is the same as in panel (a). 


can infer the presence of a finite electric field; (iii) the 
asymptotic limit of |Af| is still given by 2e/(3a). 

In the case of the S-dimer configuration, the bond be¬ 
tween S atoms leads to a reduction in the electron affin¬ 
ity of the dimer with respect to atomic sulfur. As a 
consequence, the dimer is able to only partially trans¬ 
fer the electrons that for bare edges would be captured 
by the outer Mo atom. This means that the +2e bound 
charge per unit length associated with the empty Wan- 
nier function is no longer fully screened by the nega¬ 
tive charge taken by the S dimer. Thus, one has that 
A^° = Ap° + 2e/a = 8e/(3a), while at the S-edge 
~ ~ —2e/(3a). An electric field due to these 
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EIG. 9: Top view of a 112-atom M 0 S 2 triangular island with 
zigzag Mo-edges saturated by S dimers. As in Eig. 1, Mo and 
S atoms are represented as dark grey and light yellow cir¬ 
cles, respectively, while small red circles denote the positions 
of the Wannier functions centers associated with the occu¬ 
pied electronic states. Inset: Schematic representation of the 
triangular island with vectors representing the directions of 
the three possible polarization vectors with smallest magni¬ 
tude which are compatible with the bulk polarization lattice 
in Eq. (2). While in the interior (white region) all vectors are 
equivalent, close to the edges (colored regions) the presence of 
a specific termination constrains the polarization to a single 
physical value. 

bound charges is still present, as shown in Fig. 8(b), and 
points in the same direction as for the case of bare edges. 
We note that at the Mo-edge the spatial separation be¬ 
tween the partially screened polarization charge and the 
negative charge accumulated on the S dimer gives rise 
to a finite dipole signaled by a locally negative slope in 
v{x). A charge reconstruction occurs to screen the bound 
charges and the system becomes metallic (see panel (a) 
in Fig. 8). In agreement with Ref. 5, three bands cross 
the Fermi energy: one is associated with the metallic 
states on the S-edge and it is unaffected by the presence 
of the dimers; the other two are related to metallic states 
on the Mo-edge, one having contributions mainly from 
the S dimer while the other from outmost Mo atoms. In 
this case the asymptotic limit of the free charge density 
is different at the two edges since the bound charge is 
different, i.e. |A^°| approaches 8e/(3a) asymptotically, 
while |Ap| ^ 2e/(3a). 

What has been discussed so far for zigzag nanoribbons 
is also key to understand the presence of edge states 
also in triangular islands^’^^. We show in Fig. 9 the 
relaxed structure of a 112-atom M 0 S 2 triangular island 
with zigzag Mo-edges saturated by S dimers. We will 
focus on this edge structure since according to previous 
calculations it is the most stable at standard experimen¬ 
tal conditions^’^^, even though the following arguments 
apply also to different edge terminations. Since we are 
now dealing with a finite-dimensional system we can de- 


















fine the total electric dipole, which can be computed e.g. 
by mapping the occupied states into maximally-localized 
Wannier functions as was done in Eq. (1). The centers of 
these Wannier functions centers are reported in Fig. 9 as 
small red circles. In agreement with three-fold rotation 
symmetry, the total electric dipole vanishes. Nonethe¬ 
less, this does not mean that locally the electric dipole is 
identically zero. In addition, we note that in the interior 
of the triangular island the Wannier functions’s centers 
are in the same positions as in bulk M 0 S 2 (see Fig. 1). 
From the earlier discussion, we recall that the bulk po¬ 
larization of group-VI TMDs is always non-zero, and is 
given by a lattice of vectors generated by letting R range 
over all direct lattice vectors in Eq. (1). In particular, 
this lattice includes three possible polarization vectors 
with the smallest magnitude, whose directions with re¬ 
spect to the island edges are represented with vectors of 
different colors in the inset of Fig. 9. While in the bulk 
crystal or in the interior part of the island all these po¬ 
larization vectors are equivalent, the presence of an edge 
breaks the translational invariance and thus the equiv¬ 
alence between the different polarizations (see Ref. 44 
for a related discussion in the case of boron-nitride nan¬ 
otubes). Indeed, a specific edge termination reduces the 
freedom in the choice of the bulk unit cell and thus con¬ 
strains the possible values of the polarization which are 
consistent with the crystal termination^^This was 
implicitly assumed in Eq. (2) where a specific value of 
the polarization was adopted according to the nanorib¬ 
bon termination. In the case of triangular islands we have 
three different edge orientations, each compatible with a 
different polarization vector which is orthogonal to the 
edge. We can thus say that, while in the interior of the 
triangular island all polarization vectors are admissible, 
close to the edges only a specific polarization emerges 
as schematically represented in the inset of Fig. 9. This 
preserves the zero total electric dipole but gives rise to 
a finite (positive) polarization charge at the edge of a 
triangular island. A charge reconstruction then occurs 
in order to screen this polarization charge (together with 
possible additional bound charges) with the accumula¬ 
tion of electrons close to the edge and the emergence of 
mid-gap states close to the chemical potential. Thus, also 
in the case of triangular islands, we can explain the exis¬ 
tence of such mid-gap edge states (already identified both 
theoretically^ and experimentally^^) as a consequence of 
the polar discontinuity occurring at its boundaries. A 
connection with the case of nanoribbons can be drawn in 
the asymptotic limit of infinitely large triangular islands, 
when each edge is independent of the others and the sys¬ 
tem becomes metallic with a continuous set of mid-gap 
states localized close to the edges and with a charge equal 
and opposite to the bound charge. 

Zigzag nanoribbons and triangular islands are not the 
only expression of metallicity in group-VI TMDs. Indeed, 
as we mentioned above, there is both theoretical^^’^^ and 
experimental^^ evidence that metallic states appear also 
at zigzag inversion domain boundaries (IDEs). These 


grain boundaries emerge as a result of the lack of in¬ 
version symmetry in group-VI TMD monolayers when 
two crystallites related by inversion are merged together. 
An explanation of the origin of metallicity in zigzag IDE 
is still missing. In the following we show that when 
the IDE extends along a zigzag direction a polar dis¬ 
continuity emerges across the boundary where polariza¬ 
tion charges appear. Eoundary states are then induced 
to screen these polarization charges and the system be¬ 
comes metallic. In order to verify this picture within first- 
principles DFT simulations we considered supercells con¬ 
taining two IDEs extending along a zigzag direction and 
separating two stripes of M 0 S 2 with opposite atomic ar¬ 
rangement. We assumed the structure of the S-boundary 
to be the one identified in experiments^^ (see the right 
boundary in Fig. 10(a)). In order to assess the struc¬ 
ture of the Mo-boundary we performed several calcula¬ 
tions for different atomic arrangements, assuming both 
stoichiometric and non-stoichiometric configurations. In 
agreement with Ref. 21, we obtained that the stoichio¬ 
metric Mo-boundary structure shown in Fig. 10(a) is the 
most stable over a wide range of sulfur chemical poten¬ 
tials and we will assume it in the following. We note that 
for this structure both the S- and Mo-boundaries display 
reflection symmetry through a vertical plane containing 
the interface. The energy bands as a function of momen¬ 
tum along the IDEs are plotted in Fig. 10(b) and are 
consistent with the results reported in Ref. 21. In partic¬ 
ular, we note that the system is metallic with two bands 
crossing the Fermi energy: one is a valence band that 
gets partially depleted while the other is a partially-filled 
conduction band. The charge density of free carriers as¬ 
sociated with these two bands is shown in Fig. 10(a), 
with holes localized on one IDE and electrons on the 
opposite IDE. Metallicity is induced by the electrostatic 
potential -whose macroscopic average v{x) is shown in 
Fig. 10(c)- that shifts the bands upwards in energy at 
one IDE creating a hole pocket, while downwards at the 
opposite IDE giving rise to an electron pocket. To un¬ 
derstand the origin of this electrostatic potential we first 
recall that because of the opposite orientation of the crys¬ 
tals facing the IDE, we have a polar discontinuity across 
the boundary. For an IDE extending along a zigzag di¬ 
rection, we have that the density of polarization charges 
appearing at the boundary is given by Ap = ±4e/(3a), 
where the positive (negative) sign holds for the Mo(S)- 
boundary. The signs of these charges are not consistent 
with the change in slope of v{x) in Fig. 10(c), suggest¬ 
ing that there must be additional contributions to the 
bound charge at the boundaries. We first note that at 
the S-boundary, the vertical plane of S atoms is shared 
between the crystals at the two sides so that the calcula¬ 
tion of the polarization charges based on the bulk prop¬ 
erties leads to a double counting of the contribution from 
the S atoms (2 x 2e/a). In addition, the formation of a 
Mo-Mo bond leads to the localization of two additional 
electrons at the boundary (—2e/a). This means that the 
bound charge can not be identified with the polarization 
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FIG. 10: Panel (a): Top and lateral views of a M 0 S 2 structure containing two inversion domain boundaries extending along a 
zigzag direction. Green and orange triangles emphasize the opposite atomic arrangement at each side of the inversion domain 
boundaries. The real-space distribution of the charge density of free electrons (blue) and holes (red) is superimposed. Panel 
(b): Energy band structure as a function of the crystal momentum along the interface for the inversion domain boundary shown 
in panel (a).The dashed line highlights the position of the Fermi energy. Panel (c): Macroscopic average of the electrostatic 
potential energy defined in Eq. (4) as a function of the coordinate across the system. Thin solid lines identify the position of 
the two inversion domain boundaries in the supercell. Panel (d): Gharge density of free carriers |Af| as a function of the total 
width. The dashed line highlights the asymptotic limit 2e/(3a). The solid line is a fit according to Ref. 24. 


charge reported above but has to be corrected in the fol¬ 
lowing way: = Ap + 4e/a — 2e/a = +2e/(3a). At the 

Mo-boundary, we have instead an extra row of S atoms 
which are bonded to the terminal Mo atoms at both sides. 
Owing to the different S-coordination with respect to the 
bulk, six electrons are needed to form such bonds. Four of 
them come from a charge reorganization at the interface 
and would otherwise be associated with Wannier func¬ 
tions localized at the boundary. The last two supplemen¬ 
tary electrons give rise to an additional — 2e/a contribu¬ 
tion to the bound charge: A^° = Ap° —2e/a = —2e/(3a). 
We thus have that the bound charge has the same magni¬ 
tude at both boundaries and changes sign with respect to 
the bare polarization contribution. This is now consistent 
with the change in slope of v{x) and explains the origin of 
the electrostatic potential. The free carriers appearing at 
the boundaries have to asymptotically screen the bound 
charge in order to prevent a divergency in the electro¬ 
static energy. In Fig. 10(d) we plot the free charge |Af| 
as a function of the total width of the system (approx¬ 
imately twice the separation between successive IDB). 
Indeed, the asymptotic limit is consistent with the value 
of the bound charge obtained above, |Af| ^ 2e/(3a). 
We have thus identified a consistent picture to rational¬ 
ize the presence of metallic states previously reported at 
IDBs21-23. 

Before concluding we stress that these ID wires of free 
carries in group-VI TMD nanostructures have the poten¬ 
tial to deliver extremely promising applications. In ad¬ 
dition to their established role in catalysis and their rel¬ 
evance as ID channels in ultra-thin electronics and spin- 
tronics, we envision fruitful applications in solar-energy 


harvesting. Indeed, these systems support an intrinsic 
photovoltaic effect: an electron-hole pair created in the 
“bulk” semiconducting interior of such TMD nanostruc¬ 
tures is split by the built-in electric field due to the polar¬ 
ization charges. The electron and the hole are then driven 
towards opposite interfaces where they can be collected 
at the ID wires of free carries, as already discussed in 
Ref. 25. 

In conclusion, we have shown with extensive first- 
principles DFT simulations that metallicity in TMD 
nanostructures is associated with a widespread and uni¬ 
versal phenomenon driven by polar discontinuities emerg¬ 
ing across zigzag edges or inversion domain boundaries 
and it is robust over different specific edge/boundary con¬ 
figurations. The charge density of free carriers increases 
with the size of the system until asymptotically it per¬ 
fectly screens the bound charges at the edges/boundaries 
in order to prevent the divergence of the electrostatic 
potential. The polarization charges associated with the 
polar discontinuity give the major contribution to the 
bound charge at each edge/boundary. Nonetheless, the 
knowledge of the bulk polarization alone is not sufficient 
to predict the asymptotic value of the free charge density 
as additional bound charges can appear depending on the 
edge/boundary structure and induce more free carriers. 
This happens for instance when hydrogen atoms or sulfur 
dimers are chemisorbed to terminate the dangling bonds 
at nanoribbon edges. Even for finite-dimensional nanos¬ 
tructures like triangular islands, where the total electric 
dipole is zero, a polar discontinuity still occurs across 
their edges and explains the emergence of mid-gap edge 
states previously discussed in the literature. 
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Methods 

We performed first-principles density-functional- 
theory simulations using the PWscf code of the 
Quantum-ESPRESSO distribution^^. The exchange 
and correlation functional is that of Perdew-Burke- 
Ernzerhof^^ form of generalized gradient correc¬ 
tions. We use ultrasoft pseudopotentials^^ from 
PSlibrary. 0.2. 5^^ with energy cutoffs of 60 and 600 
Ry for wave functions and density, respectively. Bulk 
relaxed structures are obtained within the Broyden- 
Eletcher-Goldfarb-Shanno algorithm by requiring that 
the forces acting on atoms are below 1 meV/A and the 
residual stress on the cell is less than 0.5 kbar. Nanorib¬ 
bons are then obtained by breaking the periodicity along 
a given direction without further relaxation; only the 
position of H or S atoms saturating the edge bonds 
are relaxed in this case. We checked that including the 
atomic relaxations for all other atoms in the nanorib¬ 
bons leads only to marginal changes in the charge 
density of free carriers; only for extremely narrow ones 
{N^ < 3 for M 0 S 2 ) structural reconstructions can occur 
as a competing mechanism to compensate polarity^^. 
Einally, for inversion domain boundaries we relaxed the 
positions of atoms in the first rectangular unit cells 
close to the boundary while keeping the bulk periodicity 


along the interface. Eor nanoribbons and inversion 
domain boundaries a fine 1 x 18 x 1 Monkhorst-Pack 
grid coupled to a small cold smearing^^ (0.001 Ry) are 
used to sample the Brillouin zone and to obtain accurate 
results for the free carrier charge density. In order to get 
rid of the spurious interactions with periodic replicas 
along the vertical direction, we correct the electrostatic 
potential using the Otani-Sugino approach^^ with 20 A 
of vacuum. Eor nanoribbons, spurious periodic replicas 
are present also in-plane. In this case their electrostatic 
interaction is reduced by allowing for an amount of 
vacuum equal to the nanoribbon width (> 20 A). The 
charge density of free carriers Ap is the amount of charge 
per unit length by which valence/conduction bands 
are depleted/filled and is computed by integrating the 
corresponding density of states. 
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